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Using 226 million BB events recorded on the Y(iS) resonance with the BABAR detector at the 
SLAC e + e~ PEP-II storage rings, we reconstruct B~ — * D*°e~V e decays using the decay chain 
D*° — » D°tt° and D° — > K~tt + . From the dependence of their differential rate on w, the dot 
product of the four- velocities of B~ and D*° , and using the form factor description by Caprini et 
al. with the parameters F(l) and p\ 1 , we obtain the results p 2 Al = 1.16 ± 0.06 ± 0.08, F(l) • |V c f,| = 
(35.9 ± 0.6 ± 1.4) • 10~ 3 , and B (B~ -> D*° e -V e ) = (5.56 ± 0.08 ± 0.41)%. 
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PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



The Standard Model of particle physics (SM) con- 
tains a large number of free parameters which can only 
be determined by experiment. Precision measurements 
of all of these parameters are essential for probing the 
validity range of the model by comparing many other 
precision measurements with SM calculations. One of 
the SM parameters, the clement \V c b\ of the Cabibbo- 
Kobayashi-Maskawa quark-mixing matrix, is determined 
with semilcptonic S-meson decays. Their rates T are 
given by the universality of the weak interaction (the 
Fermi constant Gf), by quark mixing (T cx G F |V c b| 2 ), 
and by strong-interaction corrections calculated in heavy- 
quark effective QCD. For the exclusive decays B° — > 
D* + l~vi and B~ — > D*°£~V e (I = e,p), these correc- 
tions are expressed as form factors in the differential rate 
dr/dui, where w is the dot product of the four veloci- 
ties of the B and the D*. The form factors depend on 
the three parameters p 2 , i?i(l), and i?2(l) 0- Whereas 
the B° mode has been measured by many experiments 
the B~ mode has only been measured by two groups 
with much smaller data samples. However, the B° 
experiments do not agree well in their p 2 results. Using 
the isospin symmetry dT(B~ -> D* Q £-V e ) = dT(B° -> 
D* + £~v : f) , a precision measurement of the B~ mode can 
improve knowledge of p 2 and consequently of V and \V c b\- 
The aim of our analysis [H is the determination of 
the differential decay fraction dB(B~ — > D*®e~V e )/dw, 
where B = Ft, with the B~ lifetime r. The neu- 



Ha 1 (w) 



= 1 - 



trino in the B 



D 



decay is not reconstructed. 



Therefore, the w value of each reconstructed event can- 
not be obtained, only an approximation w as defined 
below. Instead of unfolding dB/dw, the parametrized 
dB/dw expectation convolved with the if resolution from 
Monte Carlo (MC) simulation is fitted to the observed 
dB/dw distribution. The fit uses the parametrization 
of Caprini et al. with p 2 = p 2 Ai and determines 
the two parameters F(l) ■ \V c b\ and p 2 . The decay 
fraction B is obtained by integrating dB/dw. Using 
the notations AAI = tub — mp* , r = mo* /tub, and 
z = (^/w + 1 — \/2) / (y/w + 1 + \/2), the parametrization 
is defined by the following expressions: 

dr Gl\V c 



dw 



F\ v cb\ 

48tt 3 



{AM) 2 m 3 D , yjw 2 - l(w> + l) 2 g(w)F 2 (w), 



g{w) = 1 + 



Aw m B — 2wmBiTijy + m D , 



w + l 



(AM) 2 



F 2 (w) = %g 1 J2 \M 



if j 



i=0,+,- 



H±(w) 



H (w) = l + ^-^[l-R2(w)}, 
1 — r 

Vl — 2wr + r 2 



1 



l z + (53p 2 - 15) z 2 - [2Z\p 2 - 91) 



Ri(w) = i?i(l) - 0.12O - 1) + 0.05(w - l) 2 , 
R 2 {w) = R 2 {\) + 0.11(w - 1) - 0.06(u; - l) 2 , 

with -F'(l) = ft-AiCl)- The values of -Ri^l) are not deter- 
mined in this analysis; they are taken from Ref. Q. 

For our analysis, we use 205 fb" 1 of e + e _ annihilation 
data recorded at y/s f=a m(T(4S)) with the BABAR detec- 
tor at the SLAC PEP-II storage rings . In addition 
to these on-peak data, we also use 16 fb _1 of off-peak data 
collected 40MeV below the T(4S) resonance. We select 
B~ — > D*°e~i/ e candidates Q by pairing electrons with 
p* > 1.2GeV/c in the e + e~ rest frame (cms) with D*° 
candidates. Since the precision of our results is not sta- 
tistically limited, we restrict the analysis to the sequen- 
tial decay modes D° — * K~ir + , which has the smallest 
combinatorial background, and D*° — > D°tt°, which has 
a better resolution in Am = m(_R'~7r + 7r ) — m(A" _ 7r + ) 
than D* a -> L>°7. 

Charged particles are selected if they have at least 
10 hits in the drift chamber, transverse momentum 
Pt > O.lGcV/c, and a polar angle between 23.5° and 
145.5° in the laboratory frame. Electrons (kaons) are 
selected with tight (loose) particle identification criteria 
[irj| . Neutral pions are reconstructed from two photons, 
each with energy above 30 MeV and a photon-compatible 
lateral shower shape in the calorimeter. The two pho- 
tons must be consistent with the ir° hypothesis (115 < 
m 77 < 150 MeV/ c 2 ). A kinematic fit with the constraint 
m 77 = m^o improves the Am resolution by a factor of 3. 
The decay candidates have to fulfill the following addi- 
tional requirements: the D*°-D° mass difference and the 
D°-candidate mass must satisfy 135 < Am < 153 McV/c 2 
and 1.8496 < m(if"7r+) < 1.8796 GcV/c 2 , respectively. 
To reject non-B-decay candidates, the second normalized 
Fox- Wolfram moment [ll[ of the event has to be smaller 
than 0.45. To help reject combinatorial background with 
a D*° and an e~ from different B mesons in the event, 
the cms angle between them must be larger than 90°. 

Since there are many low-energy background photons, 
the selection criteria result in many events with two or 
more D*°e candidates, on average 1.75 per event. All 
D*°e candidates in the same eKir combination form one 
group, called a candidate group. On average there are 
1.015 candidate groups per event. When an event has 
more than one candidate group, we keep only the one 
with the best \m{K-K) — m(D°)\. All candidates in one 
group are kept in the analysis because the simulation 
of low-energy photons is not perfect. This procedure en- 
sures that correctly reconstructed candidates are selected 
with the same probability in data and MC simulation. 

The surviving candidates are binned in Am, cos9^ Y , 
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and w. The first two variables are used for signal- 
background separation, and the third is used for the w 
dependence of the signal. The mass difference Am is de- 
fined above, and 6q Y is the angle between the B meson 



TABLE I: Summary of input parameter values. 



Input Parameter 



Value 



Ref. 



Pi 



e system in the cms defined by 
4-2(i?*i? Y -|p*||p Y |cos^ Y ) 



and the Y = D*° 
= m| + 
The value of 

w = w(/3*) = - | | cos /r)/(m B m D .) 

cannot be determined since the angle (3* between the B 
and the D*° in the cms is unknown. However, (3* is 
bounded by a minimum and a maximum value and we 
use w = [w((3* iin ) + w(f3^ nax )]/2 as an estimator for w. 
Both w and w range from 1.0 to 1.5, and the distribution 
of w — w is nearly Gaussian with an RMS of 0.026. 

The fit for V = F(l)\V c b\ and p 2 is a binned maximum- 
likelihood fit with 41, 14, and 10 equidistant bins in Am, 
cos0g Y ' an d u>, respectively. The fit function in each 
w bin is the sum of the signal function Sn,(V, p 2 ) and 
23 background functions Bi^(V, p 2 ). Each summand 
is taken as the product of one-dimensional functions of 
Am and cos0g Y . The Am distributions of correctly 
(wrongly) reconstructed D*° mesons are parametrized by 
the sum of 3 bifurcated Gaussians (product of an expo- 
nential and a power law function) . The cos 9q Y distribu- 
tions are modeled by modified KEYS functions [f|. 

The factor functions of are obtained from fits to 
the reweighted signal MC distributions with V-, p 2 -, 
Ri(l)-, and i?2( ^-dependent weights on the generator 
level. Sw also includes the total number of produced BB 
pairs, all decay fractions of sequential decays, the B~ 
lifetime, all MC reconstruction efficiencies, and efficiency 
corrections. The corrections for track reconstruction and 
charged-particle identification are obtained from control 
data samples and their MC expectations. The correc- 
tion of the 7T° reconstruction efficiency is described be- 
low. Small corrections are also applied for deviations of 
the shapes of the Am distributions in data and MC be- 
cause of track resolution differences, and for deviations in 
the shapes of the cos 6*g Y distributions because of differ- 
ences in storage-ring energy calibration and resolution. 

The background functions are separately determined 
for the 23 background classes [H]. The large number of 
backgrounds is necessary in order to factorize all B\n, as 
Bi j i ,j,(Am) x i?2.i,i«(cos 6*g Y ). The one-dimensional fit 
functions Bj\^ are again obtained from fits to MC dis- 
tributions. The fit to the data has 49 free parameters; V, 
p 2 , and 47 for adjustments of Am shapes, cos 6*g Y shapes, 
and background fractions. The number of e + e~ — > cc 
background events is fixed by the off-peak data. 

As validation of the fit procedure, we perform our fit 
on five different MC subsamples whose size corresponds 
to that of the data sample. All five results for V and p 2 
agree with the MC input to within one standard devia- 
tion. Applied to the data and using the input-parameter 



B (T(4S r ) -> B+B~ 
' -» D°7T a ) 

* K-TY+) 
77) 



B\ 




B\ 


} D° - 


B\ 


y- 






fll(l) 


fla(l) 



(50.6 ± 0.8)% [12] 

(61.9 ±2.9)% [12] 

(3.80 ± 0.07)% [12] 

(98.798 ± 0.032)% [12] 

(1.638 ± 0.011) ps [12] 

1.429 ± 0.075 [6] 

0.827 ± 0.044 [6] 



values in Table U the fit result is V = (35.9 ± 0.6) • 10~ 3 
and p 2 = 1.16 ± 0.06 with a correlation coefficient of 
+0.90. The result leads to B = (5.56 ± 0.08)% after in- 
tegrating dB/dw. The total number of signal events is 
23 499±329. A control value of x 2 can be calculated after 
the fit as a goodness-of-fit measure. We find 4436.3 for 
4095 degrees of freedom after rebinning in regions with 
low statistics. The values of % 2 in the MC-subsamplc 
fits are of similar size indicating that the factorization 
assumptions for Sw and Bi^ are not perfect. Since there 
is no bias in V or p 2 in the MC-subsample fits and no sig- 
nificant correlation between background parameters and 
both V and p 2 in the fit to the data, we conclude that 
the results are unbiased. 

Figure [1] shows the result of the fit together with the 
selected data. The "Signal" part of the fit function 
contains the correctly reconstructed B~ — * D* Q e~V e 
decays. The two D** parts contain B — > D**ev de- 
cays with ("Am peaking") and without ("Am flat") a 
correctly reconstructed D*° intermediate state (D** = 
Di, Dq, D'i, L>2, D*ir, Dw). Events with a correctly re- 
constructed D*° and a correctly identified electron from 
the same B and from two different B mesons are in the 
"Correlated" and "Uncorrelated" background parts, re- 
spectively. "Signal-like" are true decays B~ — > D*°e~T7 e 
and B° — > D* + er~D c which arc not correctly recon- 
structed. The background from true B — ► D°ev decays 
is called "D°ei/\ All other background candidates from 
BB events ("Combinatorial D*°") are flat in the Am and 
the cos #g Y distributions since they do not contain a cor- 
rectly reconstructed D*° and they do not come from a 
charmed semileptonic decay. The last contribution, only 
visible at high w, comes from cc events. 

To determine the systematic uncertainties listed in Ta- 
ble [TT] we either rerun the fit with varied input or we 
rescale the fit result. The upper part of the Table gives 
the "internal" uncertainties which are specific to our 
analysis. The relative uncertainty on the efficiency to 
reconstruct a track is 0.8%, leading to 2.4% and 1.2% for 
B and V. The dependence of the tracking efficiency on 
the transverse momentum pt has an uncertainty which 
could distort the shape of the w spectrum. The uncer- 
tainties arising from the identification (ID) of charged 
tracks as electrons or as kaons contribute to the result as 
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TABLE II: Relative systematic uncertainties in percent. 

AV/V Ap 2 /p 2 AB/B 



□ Signal 

D** (A?n-peaking) 

■ D** (Am-flat) 
Correlated 

■ Uncorrclated 



1.5 
w 

Signal-like 
D°eu 

Combinatorial D*° 
cc events 



FIG. 1: (Color online). Data distributions (dots with er- 
ror bars) and fit results (stacked histograms) for (a) Am in 
the cos^by signal range (— 1,+1), (b) cos#b Y m the Am sig- 
nal range (140,144 MeV/c 2 ), and (c) w in both signal ranges. 
The plot below (c) shows the quotient fit/data. The different 
contributions to the fit function are explained in the text. 



listed under "particle ID efficiency". A significant frac- 
tion of the total uncertainty comes from the precision of 
the 7T° reconstruction efficiency (e^ )- It is determined 
from e + e~ — > t + t~ events where one of the two r lep- 
tons is either reconstructed by one track and two clusters 
(mainly r — > p^tt )^) or by only one track without clus- 
ters (mainly r — > iris, p,uV), The other r, used as a T-pair 
tag, is reconstructed in its evv decay. From the numbers 



Tracking efficiency (etr) 


1.2 




2.4 


TJti (\ PnPTlHpTIPP Pif" f+ 


0.3 


0.5 


0.2 


Particle ID efficiency 


0.9 


2.0 


1.6 


Extrapolated it efficiencv (e„o) 


1.8 




3.6 


p n o dependence of t n o 


1.0 


3.5 


0.4 


Am shape of D** background 


0.1 


0.1 


0.2 


Shape parameters 


1.0 


2.5 


0.6 


Number of 13 H events 


0.6 




1.1 


Off-peak luminosity 


0.1 


0.4 


<0.1 


MC statistics 


0.3 


0.8 


0.2 


Radiative corrections 


0.5 


0.4 


1.4 


Total internal 


2.9 


4.9 


5.0 


Ri(l) and ife(l) 


0.4 


4.7 


0.5 


B(T(4S) -> B+B~) 


0.8 




1.6 


B (D*° -» D°7T ) 


2.3 




4.7 


B(D° -> A'~7r+) 


0.9 




1.8 


_B~ life time 


0.3 






_D** decay fractions 


0.3 


0.7 


0.3 


Number of D*° in cc events 


0.2 


0.7 


<0.1 


Total external 


2.7 


4.8 


5.3 


Total 


3.9 


6.8 


7.3 



of t + t~ events reconstructed in each of the two chan- 
nels we derive an efficiency in data and in MC, giving 
a correction to the simulated ir° efficiency. The correc- 
tion is obtained for momenta above 350 McV/ c and has a 
precision of 3%. In the lower- momentum region with all 
7T° mesons from D*°eis decays, we use a correction factor 
of 0.960 ± 0.035 where the increased uncertainty covers 
the extrapolation into this region. Efficiency differences 
between t + t~ and BB events are covered by the MC 
simulation as controlled by comparing the rates of recon- 
structed D° decays into K~ ir + and K~ tt + it . The 
uncertainty in the shape of the w spectrum, i.e. its influ- 
ence on p 2 , is estimated by fit results for different lower 
cutsonp^o ("p T o dependence of e T o"). Corrections to the 
Am shape and to the cos 9^ Y shape are parametrized as 
functions of w, see "shape parameters" for their contri- 
butions to the systematics. Uncertainty estimates from 
radiative corrections are taken from the BABAR analysis 
of B° — > D*lv decays @ which uses the same lepton- 
momentum cutoff of 1.2 GeV/c. 

The "external" uncertainties owing to parameters 
taken from other experiments arc given in the lower part 
of Table [ill For p 2 they are dominated by Ri(l) and 
i?2(l)- For future updates, we also give in Table [DTI the 
derivatives of our three results with respect to these two 
variables as determined from fits with varied input val- 
ues. The B — > D**ev decays contribute to the uncertain- 
ties because of their less precisely known decay fractions 
and their uncertain Am shape due to low-energy photon 
background. Uncertainties in their w shape are covered 
by 10 of the 49 fit parameters. 
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TABLE III: Derivatives of V, p 2 , and B. 





V 


P 2 


B 


8/dR^l) 


-0.00342 


+0.0303 


-0.00567 


d/dRaQ.) 


-0.00525 


-1.22 


-0.00594 



Adding all systematic uncertainties in quadrature leads 
to the last line in Table Hi] and to our final results 



F(l) ■ \V C 



cb 



(35.9 ± 0.6 ± 1.4) ■ 10" 3 
1.16 ±0.06 ±0.08 , 



B{B~ -> D* e-V e ) = (5.56 ± 0.08 ± 0.41)% . 

The correlation coefficients between F(l) • \V c b\ and p\ 
arc +0.90 for statistics, +0.42 for systematics, and +0.52 
in total. Using F(l) = 0.919 ± 0.033 from lattice QCD 
0, we obtain |F c6 | = (39.0 ± 0.6 ± 2.0) • 10~ 3 in good 
agreement with the average from the exclusive neutral 
B decays B° -> D*~£ + v, (39.2 ± 0.7 ± 1.4) • 10" 3 0, 
and in agreement with results from the inclusive decays 
B -> X c £u, c. g. (42.0 ±0.2 ±0.7) ■ 10" 3 in Ref. Q. Our 
result for p 2 is in the center of the range (0.5, 1.5) from 
the B° -> D*-(+v experiments 0. 

Compared with the PDG average [l2| of B(B~ — > 
Z)*°e~z7 e ), our result is lower by more than 1.5 standard 
deviations. For a comparison of our decay- fraction re- 
sult with that of the B° mode, we use t(B+)/t(B°) = 
1.076 ± 0.008 and B(B a D*-£+v) = (5.28 ± 0.18)% 
0. This gives B(B~ -> D*°£-V) = (5.68 ± 0.20)%; our 
result agrees well with this value. 

To conclude, this measurement is the first one for 
B~ — > D* £~Ve decays with a data sample comparable 
to recent B° — > D* + £~vg experiments. The results for 
the decay rate and for \V c b\ agree well with the B° mean 
values. Since the uncertainties in the reconstruction of 
low-momentum ir + and ir° are experimentally very dif- 
ferent, the agreement of our p 2 result with the central 
value of the B° results provides a crucial cross check for 
previous \ V c b\ determinations in B — > D*£i/g decays. 
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